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Abstract. We study the deviations from perfect memory in negative temperature cycle spin glass exper- 
iments. It is known that the a.c. susceptibility after the temperature is raised back to its initial value is 
superimposed to the reference isothermal curve for large enough temperature jumps AT (perfect mem- 
ory). For smaller AT, the deviation from this perfect memory has a striking non monotonous behavior: 
the 'memory anomaly' is negative for small AT's, becomes positive for intermediate AT's, before vanishing 
for still larger AT's. We show that this interesting behavior can be reproduced by simple Random Energy 
trap models. We discuss an alternative interpretation in terms of droplets and temperature chaos. 

PACS. 75.50.Lk Spin glasses and other random magnets - 05.70.Fh Phase transitions: general studies - 
64.70.Pf Glass transitions 



1 Introduction 

It is well known that in glassy systems, dynamical effects 
strongly depend on the history of the system after quench- 
ing from above the glass transition temperature Tg . These 
phenomena are called aging and have been studied using 
various experimental protocols The measurement 

of the ac-susceptibility during negative T-cycles is one 
of them. This experiment consists of the following three 
stages. In the first stage, the system is quenched from 
above its critical temperature Tg and it is kept at a tem- 
perature Ti (< Tg) during a time t^i. In the second stage 
the temperature is temporally reduced to T2 = Ti — AT 
during a time tw2 , and then it is set back to Ti in the third 
stage. The ac-susceptibility (magnetic, dielectric, mechan- 
ical,..) is measured during all the three stages. The effect of 
the perturbation of the temperature is examined by com- 
paring the perturbed and unperturbed (i.e. tw2 = 0) data 
in the third stage. From this comparison, it is revealed that 
the perturbed data quickly approaches the unperturbed 
one as if the system remembers how far the relaxation at 
the temperature Ti had proceeded before the perturba- 
tion, even though the system is strongly rejuvenated at 
temperature T2 (see Fig. |l|). This phenomenon is called 
the memory effect, and has been observed in many glassy 
systems like spin glasses |H,^,^|^, orientational glasses]^ 
polymer glasses jTj7 etc. 
In the present paper, we focus on the deviations from 
perfect memory that are observed immediately after heat- 
ing back the system to Ti . Surprisingly, systematic exper- 
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iments (that we report in section ^) show that the tran- 
sient behavior is non monotonous as a function of AT. 
The initial extra contribution, that we will call the mem- 
ory anomaly, is found to be negative for small AT, then 
positive for intermediate AT, before vanishing completely 
for large enough AT (perfect memory). 

We then measure the ac-susccptibility for a T-cycle in 
the Random E nerg y Model (REM) Q or the Generahzed 
REM (GREM)(l3|. These models have been shown[|l|,|l|, 
to reproduce many of the experimental fea- 
tures of aging, including rejuvenation and memory. We 
find that the non monotonous transients mentioned above 
can also be obtained in such models. 

Technically, we first establish a relation between the 
ac-susceptibility and the distribution of relaxation times. 
Although this relation could not be derived analytically, 
numerical tests suggest that this relation holds with high 
accuracy. This enables us to measure the ac-susceptibility 
for any desired time scale in the REM (but not in the 
GREM), allowing us to measure the ac-susceptibility in 
the REM up to time scales comparable to experiments. 
This would not be possible using Monte Carlo simula- 
tions. As for the GREM, the ac-susceptibility is obtained 
using Monte Carlo simulation, and therefore corresponds 
to rather small time scales. Our numerical measurements 
are made for various sets of parameters, i.e., the waiting 
times iwi and iw2, the temperatures Ti, T2 and the period 
of the applied ac- field P — 27r/[j. We find that the non 
monotonous effect described above depends quite sensi- 
tively on some of these parameters. 



2 



M. Sasaki et al: Deviations from Perfect Memory in Spin Glass Temperature Cycling Experiments 




0.040 ■] ■ 1 • 1 ■ , • , 

10000 20000 30000 40000 

time [s] 



Fig. 1. Out-of-phase susceptibility x" vs. time for tempera- 
ture cycling experiments at Ti = UK on the CdCr1.7Ino.3S4 
Heisenberg spin glass (Tg — 16.7-?^). The sample was quenched 
to Ti, kept at this temperature during ti — 7700s, and then 
submitted to a negative temperature cycle at Ti — AT during 
t2 = 23650s before returning to Ti for a time is. Four different 
small AT (O.IK-OAK) were used. 
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Fig. 2. Superposition of x"-relaxations of Fig.l after the neg- 
ative temperature cycles onto a reference isothermal relaxation 
curve (solid line). The data measured after the cycle have been 
shifted horizontally by t^ — taH to take into account the effective 
contribution tcff of aging at Ti — AT on aging at Ti . The merg- 
ing of the data points with the reference curve occurs first from 
below for small AT's and then from above for larger AT's. 



The organization of this manuscript is as follows: In 
section ^ results of new complete set of T-cycle experi- 
ments in spin glasses are shown. In section^ we explain the 
REM, the GREM and the dynamics employed for these 
models. In section ^ a relation between ac-susceptibility 
and distribution of relaxation times is proposed and its 
validity is tested numerically. In section |^ and ^ results on 
the REM and on the GREM are shown. Finally, in section 
, we give a physical discussion of our results and a com- 
parison with the scenario of temperature chaos, where the 
existence of an overlap length is assumed. 



2 Results of Experiments 

The effect of temperature changes on aging has already 
been largely investigated experimentally in spin glasses Q . 
Here we focus on the details of the experimental results 
in the well characterized thiospinel CdCr1.Tlno.3S4 {Tg = 
16.7i^) Heisenberg spin glass sample. Fig. |l| presents the 
results on the out-of-phase component of x" of the a.c. 
susceptibility during the temperature cycle described in 
the introduction, with Ti = 14K, for four AT values in the 
range (O.IK-OAK), and a frequency 0.1 Hz. Note that the 
period of the a.c. field (10 sec) is typically much larger 
than any microscopic time scale (10~^^ sec). |] 

Just after the initial quench, x" is slowly relaxing down- 
ward with time twi due to aging. When Ti is decreased 
to T2 = Ti — AT, we observe a jump and a strong relax- 
ation in x"'- this is the rejuvenation effect. Despite this 
strong reinitialization of aging at Ti — AT, it is possible, 

^ The relevant 'microscopic' time scale may however be 
strongly renormalized by critical fluctuations, as recently dis- 
cussed in refs. [Q^,^, and can be much larger than 10"^'^ 
sec. for T close to Tg. 



when the temperature is raised back to Ti and for large 
enough AT (> 2K) Q, to find a perfect memory of the 
past relaxation at Ti. For large Z\T's, the relaxations at 
Ti before and after the temperature cycle are in exact con- 
tinuity and there is no contribution of aging at Ti — AT 
on aging at Ti. In contrast, in the regime of small AT^s of 
Fig. |l|, we do not find a perfect memory of aging after the 
temperature cycle. The x" relaxations after the negative 
cycle can still be superposed, apart from a transient con- 
tribution to be analyzed below, onto a reference isother- 
mal relaxation at Ti but we now need to shift the data 
by an effective time <cff < iw2 which accounts for aging 
during the stay at T2 (see Fig. ||). This effective time 
has been recently studied in detail both experimentally[^ 
m and numerically HH. 

Coming back to the transient contribution, we see that 
for the smallest AT = O.IK and 0.2K used, x" reaches a 
maximum as a function of time, and merges back with 
the reference curve from below, while for AT = Q.3K 
and 0.4^ the maximum disappears and this return occurs 
from above. This is a systematic effect which is also ob- 
served at other temperatures. We have further character- 
ized this feature by gathering the results of several nega- 
tive temperature cycling experiments done on the thiospinel 
sample at two temperatures Ti = 12K and Ti = 14K and 
for various Z\T's. In Fig. |[ we have plotted the (relative) 
memory 'anomaly' Ax" /x" j where Ax" is the difference 
between the x"- value just after the cycle and the one cor- 
responding to an isothermal aging at T during t^i +ioff , as 
a function of AT/Ti for the whole set of available experi- 
mental results. The characteristic time scales are given in 
Tables 1 and 2. 

For small Z\r's, Ax" /x" is negative and x" approaches 
the reference curve from below. As AT increases, this ratio 
becomes positive meaning that the approach takes place 
now from above the reference curve (Fig. |^). For larger 
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Fig. 3. Relative difference (memory 'anomaly') Ax" /x" be- 
tween the x"-value after a negative temperature cycle of am- 
plitude AT at a temperature T and the x"-value after a cor- 
responding isothermal aging at T during a time ti + tefi{AT) 
for various temperature cycling experiments (such as the ones 
presented in Fig. 1-2). The first cycles at lAK and 12K 
(O.lHz) correspond to ti = 7700s and t2 = 23650. The oth- 
ers (*) correspond to U = t2 = 3600s for 14K (OAHz) and 
= t2 = 18000s for 12K (O.lHz). The thin line is a guide for 
the eye. The dotted line corresponds to Ax" /x" = -AT/T. 
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Fig. 4. Structure of the Random Energy Model. 

The time t* at which the maximum of x" occurs rapidly 
decreases as AT increases, whereas the recovery time tree, 
at which the signal merges with the shifted reference curve 
has a non monotonous behavior with AT (see Tables 1 and 
2). For still larger ZiT's, this time decays back to zero. In 
addition, quite surprisingly, t^cc. is very long for interme- 
diate ZiT's - much longer than expected from activated 
thermal slowing down. This second non monotonous be- 
havior of tree, and the unexpected long ti-cc. will be recon- 
sidered in section ^ 



Table 1. tcs'- shift time, t*: position of the maximum of x", 
ticc.'. time at which the signal merges with the shifted reference 
curve, and Ax" /x"'- amplitude of the memory anomaly for 
different AT's, for Ti = 14 K. 



AT (K) 


tcs (sec) 


t* (sec) 


froc. (sec) 


Ax"/x" 


0.1 


12000 


1800 


4200 


-2.4 X 10"^ 


0.2 


4000 


900 


900 


-2.6 X 10"^ 


0.3 


1500 


160 


8000 


2.8 X 10"^ 


0.4 


1000 




10350 


9.1 X 10"^ 



Table 2. Same caption as Table 1 for Ti = 12 K. 



AT (K) 


feff (sec) 


t* (sec) 


tree, (sec) 


Ax"/x" 


0.1 


15000 


1200 


2850 


-1.8 X 10"'' 


0.2 


8000 


1200 


2700 


-2.3 X 10"^ 


0.3 


4000 


800 


800 


-1.7 X 10'^ 


0.4 


2700 




5300 


2.0 X 10"^ 



Z\T's, Ax" /x" shows a maximum and decreases back to- 
wards zero. Beyond that point, rejuvenation and full mem- 
ory effects are observed and aging at T2 has no influ- 
ence on the aging at Ti. This characteristic oscillation of 
Ax" /x" as a function of AT is the central result of this 
paper, that we discuss below in the context of Random En- 
ergy Models. It is worth noticing that in the measurement 
of d.c. magnetization with T-cycle, the relaxation rate S{t) 
at a short fixed time show a similar non monotonous be- 
havior on Z\T[M. This is quite reasonable if we notice the 
rough relation x^'C'jj) S{uj^^) (we thank to Nordblad for 
telling us this point). 



3 Models 

This section is devoted to introducing the REM, the GREM 
and their dynamical extensions. A magnetization-like vari- 
able is introduced in order to define and measure an ac- 
susceptibility. 

3.1 The REM 

The REM is schematically shown in Fig. ^. The bottom 
points are the accessible states of the system. We will con- 
sider the case where N, the total number of states, is very 
large. Each branch represents the barrier energy E, over 
which the system goes from one state to another. The 
values of E are assigned randomly and independently ac- 
cording to the distribution : 

dE 

p{E)dE = — eM-E/T,] {E > 0), (1) 

where Tg is the transition temperature of the model. Here- 
after Tg is used as the unit of temperature and is set to 
1. 

From the Arrhenius law, the escape time T(a) is re- 
lated to E{a) as 

Tia) ^Toexp[Eia)/T], (2) 

where T is the temperature and tq is a microscopic time 
scale. Hereafter tq is used as the unit of time and is set to 
1. From eq. (^, the distribution of r is 

PAr)dT = ^dr (r > 1), (3) 
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Fig. 5. Structure of the Generalized Random Energy Model 
with L — 2. 



where x = T/Tg. From eq. (||), it is clear that the averaged 
relaxation time is a;/(x — 1) for a; > 1 and infinite for x <\. 
This means that the transition from an crgodic phase to 
a non-ergodic phase occurs at Tg[I^. 

We define the dynamics of the REM from a simple 
Markoff process that defines a 'trap' model (see also ||2^). 
At t = 0, an initial state (3 is chosen according to the 
uniform distribution over all states, i.e.. 



Pa(0) 



1 

TV' 



(4) 



This means that the system is quenched from an infinitely 
high temperature. After the initial state is chosen, the sys- 
tem successively changes its state by repeating the follow- 
ing two processes. 

1. The system is activated from the present state /? with 
probability t(/3)~^ per unit time. 

2. After the activation from (3, the system falls in one of 
all the states with uniform probability. 

When a magnetic field H{t) is applied, the energy of a 
state j3 is shifted by —H{t)Mp and the activation energy 
changes from E{(3) to E{(3) + H{t)Mj3, where M/s is the 
magnetization of a state /3. This is the only effect of the 
magnetic field that we consider. The values of the mag- 
netizations Mp are assigned randomly and independently 
from a given distribution D(M) with zero mean. 



Therefore, in this model, the system freezes progressively 
from the uppermost (the L-th) layer to the lowest one as 
the temperature decreases. 

Now let us turn to the dynamics of the GREM. The 
initial state is given in the same way as the REM, i.e., 
eq. (^). After the initial state is chosen, the system suc- 
cessively changes its state by repeating the following two 
processes. 

1. The system is activated from the present state f3 to its 
A:-th ancestor /3fe (see Fig. ||) with the probability 



WiP;k) 



n 



exp[-£;,(/3,_i)/T] 



l-exphSfe+i(/3fe)/T] 



(7) 



per unit time. By convention, El^i{[3l) = oo. The first 
factor in the right hand represents the probability that 
the system is activated from (3 to (3k and the second 
one insures that the transition from (3k to /3fe+i is not 
active. 

2. After the activation from (3 to /J^, the system falls to 
one of all the states "under" (3k with uniform proba- 
bility. 

When magnetic field H{t) is apphed, Ei{f3) in eq. @) 
is replaced by Ei{(3) + H{t)Mf3. In order for nearby states 
to have strongly correlated magnetizations, the value of 
magnetization of a state (3 is assigned to be 

Mf3=MiiP)+M2{(3i) + --- + ML{(3L-i), (8) 

where A4k+i{Pk) is the contribution from the branching 
point (3k- The value of A^fc is assigned independently and 
randomly from a given distribution Dk{M.k) with zero 
mean. If the distance d{a,(3) between a and (3 is fc, i.e., 

— (3n for n > k, the correlation between Ada and 
comes from the common contributions of A4„ (n > fc -I- 1) 
to these magnetizations, and is given by 



M^Mb = 



L 

E 

n=k+l 



Ml 



(9) 



3.2 The GREM 

The GREM is schematically shown in Fig. ||. This model 
is generated by piling up L different REM's in a hierar- 
chical way. The Energy of a branch in the n-th layer (n is 
counted from the bottom). En-, are given according to the 
distribution 

AE 

p„{En)dEn^ —f-e^^[~En/Tg{n)\ (£;„ > 0). (5) 

The transition temperatures for the layers are chosen so 
as to satisfy 

Tg{l)<Tg{2)<---<Tg{L). (6) 



where A^^ is the variance of I?„(A^„). It decreases mono- 
tonically as k increases and thus as the barrier between 
the two states becomes higher, just as occurs in the SK 
model 0- 

4 Estimates of the a. c. -susceptibility 

Before we show our results, let us explain how we mea- 
sure the ac-susceptibility. One simple way is to perform a 
Monte Carlo simulation. But time scales for which we can 
study by Monte Carlo simulation is quite different from 
the experimental one. For example, if we measure the pe- 
riod of the applied ac-field in units of the microscopic time 
of the system, a typical value in numerical studies is 10^ 
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Fig. 6. The REM out-of-phase ac-susceptibility after 
a negative T-cycle. It is measured in two different ways. One 
can measure x"{^,'t) directly by a Monte Carlo simulation (the 
lines). The another way is to measure using the rela- 

tion eq. (|l^) (the diamonds) . After the system is quenched from 
an infinitely high temperature, the temperature is changed as 
Ti = 0.6 ~*T2 = 0.4 Ti. The period of the applied ac-field 
P is 100. For the data obtained by Monte Carlo simulation, an 
average is taken over 2 x 10^ samples. 



(see refs. [|l6|,^|2^,g^,^), while that in experiments is 
10^ — 10"'^^ (in the vicinity of Tg, the microscopic time scale 
may be renormalized by critical fluctuations: see refs. m, 

li). ^ 

In the case of the REM, we can overcome this problem 
by using the relation 



M2 



1 



1 — iujT 



Qir,t)dT, 



(10) 



where uj is the angular frequency of the ac-field, Af ^ is the 
variance of D{M) and Q(T,t) is the probability density 
that the system is found at time t in one of the states 
whose relaxation time is r. In eq. (10), both x(w,t) and 
Q{T,t) are the disorder averaged quantities. The merit of 
this method is that we can calculate Q{T,t) at arbitrary 
time t for arbitrary initial condition Q{T,t = 0) because 
the Green function Gpa{t), i.e., the probability that the 
system which initially is at a reaches (3 at time t, has 
already been calculated analytically in ref. jl^]. As a re- 
sult, we can estimate ac-susceptibility even for very long 
time scales comparable to those in experiments. The de- 
tails of how we can calculate Q{T,t) are described in the 
appendix. 

Here the question is whether the relation eq. ( pO[ ) is 
valid or not even if the system is not equilibrated. In 
order to examine this question, we compared data ob- 
tained by Monte Carlo simulation and those obtained by 
eq. (|l^). One example is shown in Fig. ||. A negative T- 
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Fig. 7. The REM out-of-phase ac-susceptibility x" after a 
negative T-cycle, plotted as a function of t/P, where r is the 
time elapsed after the temperature is returned to Ti and P is 
the period of the applied ac-field. After the system is quenched 
from an infinitely high temperature, it is kept at Ti = 1.1 for 
tw = 10^*. Then the temperature is reduced to Ti — AT for 
tw, and is then shifted back to Ti. The period of the applied 
ac-field is P = 10^". 



cycle is applied during the measurement. The agreement 
between both data is almost perfect. We also checked 
that both data coincide very well also for the in-phase 
ac-susceptibility x' ■ We can therefore trust the validity of 
eq. (10) for our purposes. 

Concerning the GREM, xi^^t) has been measured by 
Monte Carlo simulation because we have not succeeded in 
calculating (5(t, t) analytically. 



5 Results for the REM 

In this section, the results of ac-susceptibility measure- 
ments during a T-cycle in the REM are shown. The mea- 
surement is done in the following way: After the system 
is quenched from an infinitely high temperature, the tem- 
perature is kept at Ti in the first stage. In the subsequent 
second stage the temperature is reduced to Ti — Z\T„ 
(Z\T„ — 0.025ri), and then it is returned to Ti in the 
third stage. The time t^i of the first stage and that of the 
second stage (tw2) are taken to be equal for simplicity: 
twi = iw2 = The ac-susceptibility x(u;,t) is estimated 
by (|0|) with IP = 1. 



5.1 The case Ti > Tg 

5.1.1 Results 

In Fig. 1^, out-of-phase ac-susceptibility x" for Ti = 1.1 
and iw — 10^* is plotted as a function of t/P, where t is 
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AT 

Fig. 8. Ax" measured in the REM, plotted as a frmction 
of AT (see eq. (0) for the definition of Ax")- The values 
of the temperature Ti, the period of the applied ac-field P 



and the waiting time iw are (Ti,P,t„) = (1.1, 10^°, 10^**) 

case 1, (ri,P,tw) = (1.25, 10^°, 10^**) for case 2, {Ti,P,t^ 
^ ' (1.1, 10^", 10^3) 



(l•l,lo^^ 

case 4. 



10^**) for case 3 and (Ti, P, U 



for 



for 



elapsed time in the third stage and P is the period of the 
appHed ac-field. The value of P is 10^°. We find that x" 
approaches its equiUbrium value Xcq from below for small 
AT, it approaches from above for intermediate AT, and 
x" ~ Xcq from the beginning of the third stage for large 
AT. This is exactly what is observed in experiments. In 
Fig. H, we plot 



Ax"^x"{t = 2U + P)-x"{t = t^), 



(11) 



as a function of AT by the diamonds (case 1). We can 
easily find the similarity between this curve and the ex- 
perimental one shown in Fig. |^. 

Now let us change one of the three parameters Ti, P 
and and see how these changes affect the result. First, 
Ax" when only Ti is changed from 1.1 to 1.25 is shown in 
Fig ||by the crosses (case 2). The behavior is quite differ- 
ent from that observed in the case 1 above, i.e.. Ax" > 
for all Z\r's. Next, Ax" when only P is changed from 
10^° to 10^^ is shown by the squares (case 3). We again 
find that Ax" > for all AT's. Finally, Ax" when only 
tv, is changed from 10"^^ to 10^^ is shown by the triangles 
(case 4) . Although the value of AT above which Ax" ~ 
becomes large, there are not qualitative differences in com- 
parison with case 1. Therefore, the non monotonous tran- 
sient effect disappears (i) at low frequencies and (ii) when 
the initial temperature is not close enough to Tg. 



5.1.2 Qualitative discussion 

In order to understand these surprising results, let us in- 
vestigate the time dependent energy distribution P{E, t) 
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Fig. 9. The energy distribution P[E, t) in the REM at time 
t = 2tw-l-P after a negative T-cycle (the thick line). The values 
of parameter T, P and t„ are the same as those of the case 
1 in Fig. Sl For comparison, a function proportional to Q{E) 
(see eq. (13)) and P{E,t„) a.t t — t„ are drawn by the thin 
line and the broken one, respectively. 

which is related to Q{T,t) through the relation 

P{E,t)\AE\^Q{T,t)\dT\, (12) 

and eq. (|). In Fig. |, P{E,t) a.t t ^ 2t„ + P (i.e. in the 
third part of the cycle, after one period of the a.c. field) 
is plotted for six different values of AT (the thick line). 
The parameters Ti, tw and P are the same as those of 
the case 1 in Fig. ^. For comparison, a function which is 
proportional to 



f2{E) 



ujT _ {^)exp{E/T) 
1 + (c^r)2 " 1 + (^)2 exp{2E/T) ' 



(13) 



and P{E, t) at t = are drawn by the thin line and 
the broken one, respectively It is worth noticing that 
the data of case 1 in Fig. g are obtained by integrat- 
ing P{E,t = 2iw + P)f2{Ej over E. From this figure, 
we find that P{E,t) has both a minimum and a maxi- 
mum as a function of E (or a plateau in the special case 
Ti - Z\r = 1.0). 

When the system is kept at a temperature T for a time 
t, the equilibration at T proceeds and P{E, t) becomes 
proportional to exp[Acq(r)£'] with 



(14) 
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for < E' < riog(i). Therefore, P{E,t) is equihbrated up 
to the energy: 



E2 « (r2)iog(tw), 



(15) 



in the second stage. When the temperature is returned to 
Ti, the re-equihbration at temperature Ti starts and it 
proceeds up to 



i?3«rfog(P), 



(16) 



at time t = 2t„ + P. These considerations naturally lead 
us to the approximate shape of the energy distribution: 



P{E,2t^ + P) 



'exp[Aeq(ri)E] iO< EKEs), 
exp[Aeq(r2)E] {E3<E< E2), 



(17) 



provided E^ < i?2- Accordingly, a minimum around E w 
E3 and a maximum around E ^ E2 appear for T2 < 1.0, 
and there is a plateau between E3 and E2 for T2 — 1.0. 
As for the case T2 — 0.5, the peak is erased since in this 
case E2 < E3. 

From Fig. ^, the result of the case 1 in Fig. || is un- 
derstood as follows. Because of the existence of a peak (or 
a plateau) and the normalization of P{E,t) with respect 
to E, the difference P{E,t = 2t^ + P) - P{E,t = t^) 
changes its sign at a certain point E^, which is slightly 
greater than E3. Figure ^ shows that the peak around E2 
grows as AT increases. If we take the normalization con- 
dition into account, we notice that the growth of the peak 
means a decrease of P{E, t — 2t^ + P) for E < E^,, and a 
corresponding increase of the quantity Ax'-{AT), defined 
as 



Jo 



P{E,2t^ + P)- P{E,t^) . 

(18) 



On the other hand, if AT is not very large and the location 
of the peak is not close to i?3, the second contribution 
Ax+{AT), defined as: 



/•OO 

Ax:[{AT)^ / dEQ{E) 
Je, 



PiE,2U + P)-P{E,t„) 



(19) 



cannot be very large (note that the peak of f2{E) is around 
E^). As a result. Ax" = Ax'\. — Ax"_ is negative for small 
AT (see Fig. |l^ where Ax+{AT) and Ax'L{AT) for case 
1 are plotted by the diamonds). Then, as AT increases, 
the location of the peak approaches i?3 and Ax\{AT) 
increases. Therefore, Ax" {AT) is positive for intermediate 
AT. This peak disappears when AT is large and E2 < E^. 
This is the reason why Ax" {AT) « for large AT. 

Our next interest is to understand the following two 
trends: 

(i) If either the temperature or the period of the applied 
ac-field is large enough. Ax" {AT) > for aU AT. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

AT 



0.03 

P 

<\ 0.02 



case 

case 2 ^ 
case 3^&- 
case 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

AT 

Fig. 10. The functions Ax'U^T) and Ax-{AT) in the REM 
for each of the cases of Fig. H are plotted as a function of AT. 
See eqs. and (|l|) for the definitions of Ax" and Ax+- 

(ii) The behavior of Ax" {AT) does not depend sensitively 
on t^, except that the value of Z\r above which Ax" ~ 
increases with t^. 

To understand these trends, an important point is to note 
that Ax'-L {AT) satisfies the inequality 

/•oo 

Ax'UAT) < dEQ{E)P{E,t = U) 
Jo 

« / dEX,^{T)exp[X,^{T)E]f2{E) 
Jo 

-P^^. (20) 

It is obvious from this inequality that Ax'L (AT) is a de- 
creasing function of T and P. The quantity Ax'^{AT), on 
the other hand, does not depend strongly on either T or 
P. This explains point (i) above. 

On the other hand, there is no dependence in the in- 
equality (pO|). This is the reason why behavior of Ax" {AT) 
does not depend on so much. However, the value of AT 
above which Ax" ~ increases with increasing because 
it is determined from the condition: 



Eo ~ E^ 



(21) 



The explanation for these two trends is confirmed by Fig. 10 
where Ax+{AT) and Ax'L{AT) in each of the cases of 



8 



M. Sasaki et al: Deviations from Perfect Memory in Spin Glass Temperature Cycling Experiments 



AXl, AX" 



Tg > T2 is required to observe a non monotonous memory 
anomaly. 




6 Results for the GREM 

Now let us turn our attention to the GREM. As men- 
tioned in §^ x(^i^) for the GREM has been measured 
using Monte Carlo simulation because we have not suc- 
ceeded in calculating Q(t, t) analytically. However, we do 
not measure x(a;, t) from the linear response to an ac- field 
because this procedure requires averaging over a very large 
number of samples (typically 10^ — 10® samples). Instead, 
we have estimated the ac-susceptibility from the relations: 

Xk{u:,t) = — ^ / : Qfe(Tfe,i)drfc, (23) 



T J I - \UJTk 



where 



Fig. 11. The functions Ax'+ (the diamonds) and Ax- (the 
crosses) in the REM for the case Ti = 0.8, P = 10^° and 
tw = 10^*. The definition of Ax\. and that of Ax!_ are slightly 
changed by replacing P{E,t — t^) in eqs. (^) and ( p^ with 
Pconst{E,t = toff), where Pconst(-E,t) is isothermal energy dis- 
tribution at T and tes is defined by eq. (^) . 

Fig. I are plotted as a function of AT. The function Z\x+ (AT 
does not depend much on the parameters T, P and 
compared to Ax- (AT) , which becomes rather small when 
either T or P is large enough. 



Tk = exp 



(24) 



5.2 The case Ti < Tg 

Figure [ill shows Ax'^ and Ax'L in the case Ti = 0.8, 
P = 10^° and = 10^®. Because x" decreases towards 
zero with time for T < Tg, we change slightly the def- 
inition of AVi and that of Ax'L and replace P{E, t = 
tw) in eqs. ( p^ ) and ( |l9| ) with Pconst{E,t = tcs), where 
PconstiE,t) is isothermal energy distribution at Ti, and 
tcs is estimated as 



(iw) 



T2/T1 



(22) 



It is worth noticing that the effective time in the sec- 
ond stage is estimated as (tw)"^^^"^^- It has been shown in 
ref. that this way to estimate the effective time works 
well in the REM. The result is similar to that of the case 
Ti > Tg shown in Fig. ^ in the sense that Ax'l. has a wide 
plateau and Z\x" has a peak around the value of AT at 
which eq. (21) is satisfied (in this case the value is about 
0.356). The only difference is that /\x+ now never exceeds 
Ax'L, so that Ax" is negative for all AT. 

The dependence on the different parameters was also 
investigated by changing one of the three parameters T, 
P and The ranges we examined are 0.75 < Ti < 0.9, 
10^° < P < 10^5 and 10^® < tw < 10^^, respectively. As 
a result, we found that Ax'L > Ax'^ is always satisfied 
for all AT in all the cases. Therefore, the condition Ti > 



the function Qk (t, t) is the probability density of Tk at 
time t and Xkii^-,t) is ac-susceptibility calculated from 
Mk- It is this probability density Qk{Tk,t) that we ob- 
j;ained from Monte Carlo simulation. The validity of these 
relations was confirmed numerically by comparing data 
from the direct measurement of A4 k under an ac-field and 
from eq. (p3|). 

Because we have to rely on Monte Carlo simulation, 
the time scales are rather restricted as compared to the 
REM. Therefore, we will confine ourselves to showing re- 
sults with one set of parameters. The system we have in- 
vestigated is the GREM with L = 2 {L is the number of 
layers), Tg{l) = 0.5 and T„{2) = 1.0. The disorder av- 
erage is taken over 8 x 10 samples. The period of the 
applied ac-field is 300. After the system is quenched from 
an infinitely high temperature, the temperature is kept at 
Ti — 0.85 for tv, = 10^ in the first stage. In the subse- 
quent second stage the temperature is reduced to T2 for 
twj and then it is returned to Ti in the third stage. In 
Fig. |l^, the contribution from both levels, Xi and X2 are 
plotted as a function of r, where r is elapsed time in the 
third stage. As for x'/j we again find a non monotonous 
behavior, similar to that observed in experiments (Fig. |^) 
and in the REM (Fig. g). This was expected, since for 
the first level dynamics is very similar to the single REM, 
with transitions to the higher level frozen by the fact that 

Ti < Tg{2). 

On the other hand, the memory anomaly in X2 al- 
ways negative, for all values of AT. This result is con- 
sistent with that obtained in the REM where x" always 
approaches the reference curve from below if Ti < Tg (note 
that Ti < Tg{2) in the present case). 

Finally, let us discuss what would happen in the case 
L :i> 1. There, some layers are frozen (T < Tg{n)) and oth- 
ers layers are fast (T > Tg{n')) at any given temperature 
T <Tg = Tg{L). From the study on the REM shown in 
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Fig. 12. The GREM out-of-phase ac-susceptibility Xk after 
a negative T-cycle as a function of the elapsed time r after 
the temperature is returned to Ti. This GREM has L = 2, 
Tg{l) = 0.5 and Tg{l) = 1.0. The period of the applied ac-field 
is P = 3 X 10^ . After the system is quenched from an infinitely 
high temperature, the temperature is kept at Ti = 0.85 for 
tw = 10* in the first stage. In the second stage the temperature 
is temporally reduced to 7i — AT for the same time tw, and 
before being returned to Ti in the last stage. 



we expect that the contributions to x" from the frozen lay- 
ers will always lead to a negative memory anomaly, where 
as the contribution from the 'critical' levels will lead to a 
non monotonous contribution. 



7 Discussion. Other scenarios 

We have seen that the non monotonous transient effect ob- 
served in memory experiments can be reproduced within 
simple REM trap model, provided the temperature is above 
and close enough to the critical temperature and the fre- 
quency not too low. The same mechanism is present in 
the GREM, and is governed by the dynamics around the 
'critical level', i.e. the level such that its critical temper- 
ature is close to the working temperature. As empha- 
sized in refs. |^,^, the physical interpretation of the 
different 'levels' is in terms of length scales: small scale 
dynamics corresponds to the deepest level of the tree, 
whereas large length scales correspond to the upper level 



of the trees. The observed aging dynamics always con- 
cerns those length scales (levels) around the critical tem- 
perature: larger length scales are frozen, whereas smaller 
length scales are completely equilibrated. Hence, in the 
above GREM interpretation, the important ingredient is 
that the system remains close to criticality at any tem- 
perature, but the basic ingredient is already present in 
the REM, and is related to the abrupt change of the way 
the different states are explored at Tg (see pij). 

The GREM model is a concrete implementation of the 
so called 'hierarchical' interpretation of experimental data 
to which one often opposes the 'droplet' interpreta- 
tion . As discussed in details in ref . |^ , the two inter- 
pretations are to some extent complementary if one wants 
to interpret the 'hierarchy' of phase space as a hierarchy 
of length and time scales. 

However, the droplet interpretation of the rejuvenation 
and memory experiments makes an extra assumption that 
we now discuss. The existence of an overlap length £^t be- 
tween typical configurations at T and T—AT is postulated 
[^|3^, such that for length scales larger than I at, the 
configurations at the two temperatures are completely un- 
related ('temperature chaos'). Using plausible arguments, 
one deduces that I^t should diverge as a power-law of 
AT for AT 0. Recent experimental data has given 
some credit to the existence of temperature chaos p5[ . 
After a waiting time iwi , the active length scales are such 
that Tr{£wi, Ti) ^ twij where Tr{£, T) is the typical relax- 
ation time corresponding to length £ at temperature T. 
Length scales much smaller than are fully equilibrated. 
In this picture, the scenario for rejuvenation is thus the 
following: whenever £^i < £at, the temperature change 
does not modify the achieved pattern, but only acts to 
slow down the dynamics. Conversely, when AT is such 
that £wi > £at, the system has to start rebuilding new 
correlations as if it were brought directly from high tem- 
perature (when £at ^ 0) . As shown in details in ref. , 
this does not necessarily mean that the structure grown 
at the first temperature is immediately washed away. On 
the contrary, as long as the length scales ^w2 , active at T2 , 
remain small compared to ^wi , memory can be partially or 
totally recovered. The criterion is the following: the time 
needed to erase the effect of the dynamics at T2 when the 
system is heated back, and isothermal dynamics at Ti is 
recovered, is given by: 



^rec. (^v 



(25) 



Since £^2 decreases extremely fast with decreasing tem- 
perature [^,0,|2^,^,^, troc.(^w2) decreases very rapidly 
(for a given iw2) as AT increases, and should soon become 
smaller than ll>~^, which is the smallest time for which a 
measurement of the a.c. susceptibility can be performed. 

When ircc.(iw2) > and £^^2 > (-at, on the other 
hand, one expects to see an initial spike in the a.c. sus- 
ceptibility that corresponds to the reconstruction of small 
length scale correlations at Ti. Schematically, the tem- 
perature chaos scenario therefore predicts that the mem- 
ory anomaly Ax should be zero for AT < AT*, with 
iwi — (at* , positive for larger AT, but becoming zero 
again when ircc.(iw2) becomes shorter than uj~^. 
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One can finally argue that the number of thermally 
active (equilibrium) droplets decreases slightly when the 
temperature is reduced from Ti to T2, thereby reducing 
the equilibrium a.c. susceptibility. The need to re-nucleate 
these droplets back at Ti , which also takes a time ~ tree. , 
would then explain the negative contribution to the mem- 
ory anomaly for small AT. Q This would predict that 
Ax"/Xeq - -[l/T+\T'/r\]AT for smah AT, where T is 
the temperature dependent stiffness of the droplets. The 
experimental effect, found to be stronger than —AT/TAs 
in qualitative agreement with this prediction (see Fig. H). 
We have furthermore checked that the amplitude of the 
bump in x" at < = ^wi + ^w2 is of the same order as the 
difference between the equilibrium values of Xcq at Ti and 
Ti - AT. 

However, the time ircc. beyond which the stay at T2 
is erased does not conform to the naive estimate eq. (p5|), 
since it is found to be non monotonous and much larger 
than expected. It is rather the position t* of the max- 
imum of x" that seems to obey Eq. (p5|). Note that the 
REM scenario also predicts a monotonously decreasing re- 
covery time tree, with increasing AT . We have at present 
no physical interpretation for this discrepancy. 

The experimental data appears to be consistent both 
with the above droplet/chaos interpretation and with the 
hierarchical model developed in the present paper (see also 
the discussion in ||2^,^). The present study shows ex- 
plicitly that the non monotonous memory anomaly does 
not require the existence of an overlap length. Indeed, 
we argued that the REM trap model, where this overlap 
length is absent, is also able to reproduce qualitatively 
the memory anomaly if one works around the freezing 
temperature around which 'temperature chaos' effects are 
observed |0. In the REM scenario, the positive contri- 
bution to the memory anomaly comes from an over con- 
centration of the probability weight in deep traps at T2 
as compared to the equilibrium situation at Ti (see the 
discussion in ref. |^). Physically, this positive contribu- 
tion corresponds to a freezing at T2 of small length scales 
that have to unfreeze when back at Ti, a scenario that 
was directly confirmed by the numerical simulations of 
psf where temperature chaos is absent but rejuvenation 
and memory effects are clearly observed. 



Acknow^ledgments 

We wish to thank L. Berthier, J. Hammann, O. Mar- 
tin and H. Yoshino for many interesting conversations. We 
thank P. Nordblad and H. Yoshino for very useful com- 
ments on the manuscript. M. S. acknowledges a fellowship 
from the French Ministry of research. M. S also would like 
to thank CEA Saclay for hospitality where main part of 
this work was performed. The LPTMS laboratory is an 
Unite de Recherche de I'Universite Paris XI associee au 
CNRS. 



Appendix 

This appendix is devoted to explain in detail how we 
can calculate Q{t, t). We assume that the probability Pa{t - 
0) that the system is found at a state a at time i = is 
given. For simplicity, let us first consider the case that the 
system is kept at a constant temperature T . It is easily 
found that Q{T,t) is given as 

Q{T,t)AT = Y,Ar5{T(P)-T)Gp^{t)P^{t^Q), (26) 

where Gpa{t) is the Green function, i.e., the probability 
that the system which initially is at a reaches (3 at time t. 

Now let us calculate the Green function. When the 
system which initially is at a reaches (3 at time t, there 
are the following two possibilities: 

(i) a = [3 and the system has not been activated during 
time t. 

ill) The system is activated at i' (< t) and reaches /? after 
that time. 

In the case (ii), because the new state after the activation 
is chosen randomly from all the states, the probability that 
the system reaches (3 after the activation is Pp^^{t — t'), 
where 



puni 



(27) 



and N is the number of states. Taking this fact into con- 
sideration and recalling that the system is activated from 
a with the probability r(a)~^, we obtain 



At' 



t 



T{a) 



r(a)_ 
exp 



Pr{t-t'). (28) 



The Laplace transformation of this equation leads us to 



Giia{s) = / dtexp[-si]G/3a(i) 



ST{a) -\- 1 sr(a) -|- 1 ' 



(29) 



where P™'(s) is the Laplace transformation of P^^^it). 
From this equation and the Laplace transformation of 
eq. (27), we find 



pr{s) 



st{P) + 1 



E 



ST{a) 
ST{a) + 1 



(30) 



^ Note that the negative initial contribution to Ax is like the 
Kovacs effect in polymeric glasses |3q,b2,B3|. 



The calculation of P™'(s) for small s and its inverse Laplace 
transformation have already been done in ref. [|l5| . The re- 
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suits are 



-(/?) 



Nxs''c{x){st{/3) + 1) 

-r{P) 
iVslog(s)(sr(/5) + 1) 

{x-1){s + 1)t{P) 



(x < 1), 
(x=l), 
{x > 1), 



(31) 



and 



prit) = 



I Nxs{st{i3) + 1) 

' /q dwM^-i exp[-(t - u)/r(/3)] 
iVa;c(x)r(x) 
r(/3)[l-expH/T(/3)]] 
log(t) 

(a:-l)r(/3){l-expH/T(/3)]} 



(a:< 1), 
(x = l), 



(a:> 1), 



where x = T/Tg and 



c(x) = r{x)r{i -x)^ - 



sin(7ra;) 



(32) 



(33) 



Now let us return to the calculation of Q{T,t). The 
substitution of eq. ( p8| ) into eq. ( p6| ) leads us to 

Q(r,t) -expht/T]Q(T,0) 
df 



dr' 



'(r,i-Oe-^Q(r',0),(34) 



where 



and 



Q(t, 0)dr = ^ dr5(T(a) - T)F„(t = 0), (35) 



Q"-(r,i)dr ^ ^drJ(r(a) - r)Pr'W 

= dT7Vp,(r)P-"(T,t). (36) 

The function Px{t) is defined by eq. (^. From eqs. (^^, 
(^) and (|3^), we finally obtain 



Q(t, i) -exp[-i/r]Q(r,0) 



u — t 
exp I — — I - exp 



T' - r 

u — V 



T 



(X < 1), 



Px{t)t 
log(i) JO 



dr'g(T',0) 



|1 -e - -; |e r' 

(x-f);,.(r)r rd/g(/,o) 



{l-e-T.}--^{e-^ 



Next, let us consider how we can calculate Q{t, t) when 
the temperature is changed discontinuously as 



T{t) = T, 



{U<t< U+i). 



(38) 



'}1 (^ = 1), 



(37) 



The answer is rather simple. At first, we calculate Q{t, t^i) 
with some initial distribution Q{t,0). Then, we set the 
new initial distribution to Q(t, f^i) and use eq. ( |37| ) to 
calculate Q(r, t„2)- We can calculate Q{T,t) at any t by 
repeating this procedure. 
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